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Ever since the discovery of solitons, 1 there has been a great deal of interest in the 
existence of spatially localized solutions in nonlinear systems. In addition to the solitons of 
integrable systems, spatially localized solutions have been found in systems with dissipation 
such as the complex Ginzburg-Landau equation for systems near an inverted bifurcation 2-4 
and a nonlinear phase equation. 5,6 On a somewhat different footing, we point out that 
nonlinear evolution equations which result from perturbation expansions have been found 
to be very useful in studying the dynamics of many physical systems. An example of such 
an equation is the well-known Kuramoto-Sivashinsky equation which describes systems 
such as liquid film flow 7-9 and flame front propagation. 10 Another such equation which 
results from similar derivations 11 13 is an equation which describes the fluid-air interface 
of a thin liquid film on the surface of a horizontal plate subjected to a vertical temperature 
gradient, which gives rise to the Marangoni effect. The Marangoni effect is found to be very 
important in low gravity conditions and therefore for processes of space technology. 12 The 
one-dimensional form of the equation has been found to exhibit stable localized solutions. 13 
However, the existence of stable localized solutions for the two-dimensional equation is still 
an open question. 

In this Letter we: 1) show that the two-dimensional evolution equation, describing 
the three-dimensional spatiotemporal behavior of the air-liquid interface, has a Liapunov 
functional that is strictly nonincreasing in time and; 2) exploit this fact to demonstrate the 
existence of stable axisymmetric localized solutions for a wide range of parameter values. 
We find that for given parameter values in the equation, the stable localized solutions exist 
for a wide range of amplitudes and that there is a minimum amplitude below which the 
solution is no longer stable, but damps. Further we find that the amplitude of the localized 
solutions varies nearly linearly with the logarithm of the thickness of the film surrounding 
the localized state, larger amplitudes corresponding to a thinner film. This implies that 
the maximum amplitude of the localized state which can exist for given parameter values 


will depend on the noise in the system and on irregularities in the plate, since film rupture 
will be induced by lower noise levels and by smaller irregularities for thinner films. It will 
be very interesting to see whether experiments of a thin liquid film on the underside of a 
cooled horizontal plate will exhibit such stable stalactite-like structures. 

The equation we study is 
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Here /i(x,?/,t) is the dimensionless thickness of the film; 7, s, and B are positive dimen- 
sionless parameters describing the gravity, surface tension, and the Biot number (related 
to the rate of heat transfer at the free surface), respectively, of the system; r is a rescaled 
time; and V = ( d x ,d y ). Equation (1) results from making a thin film approximation ex- 
pansion of the incompressible fluid dynamics equations and assuming the surface tension 
is large. 

This equation was derived in Refs. 11-13 and, assuming 7 > 0 , B > 0 , and s > 0 , 
describes the spatio-temporal evolution of a thin liquid film bounded by a rigid horizontal 
plate of uniform temperature from above and by its free surface from below (the Rayleigh- 
Taylor instability); and experiencing surface tractions arising from the surface tension 
dependence on the temperature (the Marangoni or the thermocapillary effect). The value 
of the parameter m is 1 when the temperature of the plate is lower than that of the ambient 
air and —1 when the plate temperature is higher. A linear stability analysis shows that for 
m “ 1 > a hlni of uniform thickness h u is stable if h u ( 1 -fi ,B/i u ) 2 < (3 B)/ (27), demonstrating 
that the Marangoni effect can stabilize the Rayleigh-Taylor instability, and unstable for 
h u ( 1 + Bh u ) 2 > (3B)/(27), provided that the domain is sufficiently large. Further, as 
previously noted, the Marangoni effect can result in localized nonlinear drop-like solutions 
in the 1-d equation. 13 For m = — 1 the Marangoni effect is destabilizing and small unstable 
perturbations are found to evolve to film rupture. 
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Similar to the one-dimensional case, equation (1) may be written in the form 


h T = V • 




( 2 ) 


where the Liapunov functional F is 

F( r) = J dxiy [p(A) + j||Vkf] , 


( 3 ) 


and the potential P is 


p{h)=~ i -h'-+ 
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In the above and below, || |j denotes the vector norm, and the domain of integration is 
assumed to be periodic. Therefore boundary terms will vanish. 

Taking dF/dr and integrating by parts gives 


dF 

dr 
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( 5 ) 


Since h > 0, we have dF/dr < 0 and thus F is strictly nonincreasing in time. We also 
note that Q = J dxdyh is a constant of the motion. Therefore, assuming that a stable 
stationary state exists, it will be given by minimizing F subject to the constraint Q is 
constant. Taking 8{F + XQ) = 0, where A is a Lagrange multiplier, we find 



dP 

~dh. 


+ A 


(6) 


for the equation satisfied by a stationary state. This equation may be solved by finding 
the steady-state solutions of 

h ' = { V ‘ h - d T h -^ < 7 > 

where the constraint Q is constant is imposed by taking 

A (t) = -~rPr f dx dy dP/dh. (8) 

J dx dy J 
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Defining R = dP/dh - (s/3)V 2 h we have dF/dt = -JdxdyR 2 + [J dx dy R] 2 / f dx dy < 
0. Therefore, the solutions of eq. (6) found via solving the dynamical equation (7) are not 
only stationary solutions of eq. (1), but are guaranteeed to be stable stationary solutions 
since the Liapunov functional F is strictly nonincreasing during the evolution of eq. (7). 
We stress that eq. (7) does not reproduce the transient behavior of the system, but is 
only a mathematical trick to find stable stationary states. If one is interested in the time- 
dependent behavior, then it is necessary to solve eq. (1), which requires a great deal more 
computer time and is more difficult to solve due to the highly nonlinear nature of the 
equation. 

By taking m — 1, and scaling with respect to time, space, and amplitude as t = 
(8 js)/(9B 3 )T, (x,y) = y/7fi(X,Y), and h = (3B)/(2j)H, eq. (1) may be written as 

+ V • [H 3 VV 2 H] = 0. (9) 

where V = ( dx , &y ) and 0 = ( 3 B 2 )/( 27 ). To write eqs. (2-8) in the scaled variables, make 
the substitutions x -» X, y -* V, h -» H, s -+ 3, 7 -» 3, Bm -» 2 , B -» 0. We look for 
stable localized solutions of eq. (9) by numerically solving eqs. (7,8) in the scaled variables. 
We used fourth-order finite differencing for the spatial derivatives, taking a 92 by 92 grid, 
and used the Euler method for the time integration (since we are only interested in the 
stationary state, the Euler method is entirely sufficient). Periodic boundary conditions 
were used. The accuracy of the solution was checked by evaluating dF/dt from eq. (5). 

Figure la shows a stable localized solution for 0 = 1 which evolved from an initial 
Gaussian. The amplitude of this state is A = 5.7874 and the background thickness of the 
film surrounding the state is H b = 6.2789 x 10~ 2 . From eq. (8) we find A = 1.9507. Figure 
lb shows a stable localized solution with amplitude different from that of fig. la for the 
same value of 0 (note the change in vertical scale). For this state we have A = 13.6065, 
H b = 6.5057 x 10 3 , and A = 4.0542. Note that the larger amplitude state has a smaller 
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background thickness. To demonstrate that stable localized solutions exist for a wide 
range of (3 and that the form of the solutions are very similar, figs. 2a and 2b show stable 
localized states for j3 — .01 and /3 = 100, respectively (again note the vertical scale). For 
fig. 2a we have A = 12.9896, H t> — 6.0675 x 10 -2 , and A = 1.8641. For fig. 2b we have 
A = 2.1863, H\> = 3.3550 x 10“ 3 , and A = 5.2410. We see that the the width of the 
localized state is rather insensitive to the value of (3. 

Some insight may be gained by plotting the effective potential P(H) + \H as a function 
of H . Figure 3 shows such a plot for the state shown in figure la. We see that there 
is a potential well which confines the flat portion of the film surrounding the localized 
state. As H increases, the potential increases and reaches a maximum and then decreases, 
approaching — oc as H — ■> oo. In contrast to the localized states in phase dynamics 5,6 
where two potential wells exist, there is only one potential well here. Considering this fact, 
it may be surprising that stable localized states exist. However, it is the gradient term in 
the Liapunov functional that is responsible for the existence and stability of the state. 

Figure 4 shows plots of A as a function of logio(-H&) for stable localized solutions for 
5 different values of f3. To produce these plots it was first necessary to determine the 
minimum amplitude (for each value of (3) for which stable localized solutions exist. This 
was done by decreasing the amplitude of a stable state by a small amount and allowing 
the system to come to equilibrium. This process was repeated until at some amplitude 
the solution was no longer stable, but damped. Performing this process for each value of 
/ 3 then determines the endpoints at the lower right of the curves. To the lower right of 
these endpoints no stable localized solutions exist. The next step is to gradually increase 
the amplitude, beginning at the right hand endpoints and allowing the system to come to 
equilibrium at each step. As Hi became very small, it was necessary to decrease the time 
step in order to get an accurate solution (determined by evaluating dF/dt ). Therefore, 
the upper left endpoints of the curves only reflect a finite amount of computer time and 


6 



not the nonexistence of solutions. We expect that the curves continue to the upper left at 
approximately constant slope. 

An important point that can be inferred from fig. 4 is that even though stable localized 
solutions may exist for the deterministic eq. (1) for all amplitudes greater than some 
minimum, in practice - at some sufficiently great amplitude - the film surrounding the 
localized state will become sufficiently small so that rupture will be induced by noise in the 
system or irregularities in the plate. This implies that there is a maximum amplitude, the 
value of which will be determined by the noise level or by the magnitude of the irregularities. 
More precisely, we would expect that the maximum amplitude would vary approximately 
linearly with the logarithm of the magnitude of the noise or irregularities. Even if the 
plate is perfectly smooth and all external vibrations are eliminated, noise-induced rupture 
would still result from thermal noise, so that in any physical system there will always be 
a maximum amplitude for which stable localized solutions can exist. 

In this Letter we have shown that the two-dimensional evolution equation, describing 
the three-dimensional spatiotemporal evolution of the air-liquid interface of a thin film 
on the underside of a cooled horizontal plate, has a Liapunov functional that is strictly 
nonincreasing in time. We use this fact to demonstrate the existence of stable axisymmetric 
localized patterns ior a wide range of parameter values. We have found that for a given 
set of physical parameters of the system there exists a continuum of possible localized 
structures for the amplitude larger than some critical value. The initial conditions that 
can lead to such states and the transient behavior prior to the formation of the states are 
still open questions and can only be resolved upon the efficient numerical solution of the 
full two-dimensional nonlinear evolution equation. 

A.O. acknowledges the hospitality of the Institute for Computational Mechanics in 
Propulsion (ICOMP) at NASA Lewis Research Center in Cleveland where the present 



work was done. His work is partially supported by the Israel-Mexico Energy Research 
Fund and Grant-in- Aid of the Government of the State of Israel. 
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Figure Captions 


Fig. 1. Stable localized solutions of eq. (9) for f3 = 1. The arrow indicates the 
direction of gravity. Note that similar steady solutions of different amplitudes A are 
admitted by eq. (9) for the same value of /3. (a) A = 5.7874, H b = 6.2789 x 10” 2 , 
A = 1.9507; (b) A = 13.6065, H b = 6.5057 x 10~ 3 , A = 4.0542. 

Fig. 2. Stable localized solutions of eq. (9) for various values of f3: (a) (3 = 0.01, (b) 
/? = 100. Note the change in vertical scale. 

Fig. 3. The effective potential P(H) + A H for the localized state shown on Fig. la, 
as a function of H . The value of A = 1.9507 is calculated using eq. (8). The background 
thickness H b — 6.2789 x 10 -2 lies in the minimum of the potential well. 

Fig. 4. The amplitude A of the stable localized state as a function of the logarithim 
of the background thickness, log l0 (H b ) ) for various values of f3. Note that these curves are 
nearly linear. 
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(b) p = 100. 


(b) A = 1 3.6065, H b = 6.5057x1 O' 3 , X = 4.0542. 


Figure 1 .-Stable localized solutions of eq. (9) for p = 1 . The 
arrow indicates the direction of gravity. Note that similar 
steady solutions of different amplitudes A are admitted by 
eq. (9) for the same value of p. 


Figure 2.-Stable localized solutions of eq. (9) for various 
values of p. Note the change in vertical scale. 




P(H) + XH 



0.0 1.0 2.0 3.0 


II 

Figure 3.-The effective potential P(H) + XH for the localized 
state shown in Fig. 1 a, as a function of H. The value of 
X =1.90507 is calculated using eq. (8). The background 
thickness H b « 6.2789x10-2 lies in the minimum of the 
potential well. 



H b 

Figure 4.-The amplitude A of the stable localized state as 
a function of the logarithm of the background thickness, 
log 10 (H^, for various values of p. Note that these curves 
are nearly linear. 
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